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Abstract 

In the preceding paper, the exact solution of Stokes equation was obtained for a binary strong electrolyte 
solution in an external electric field. In the present paper, the solution is applied to calculate the Wien effect 
on deviation from the Coulombic law of conduction in high fields. One of the important aims of the present 
line of work was in removing or avoiding the divergence difficulty in calculating the electrophoretic and 
relaxation time coefficients. The present work achieves that aim by calculating on the basis of computing 
the axial velocity profiles the position of the center of the ion atmosphere as a function of the reduced 
field strength and therewith computing the electrophoretic and relaxation time coefficients for the migrating 
spherical ion atmosphere at each value of the reduced field strength. With the electrophoretic and relaxation 
time coefficients thus calculated along the trajectory of the center of the ion atmosphere with regard to the 
reduced field strength, the equivalent ionic conductance of a magnesium sulfate solution is calculated and 
43 . compared with experimental data in very good accuracy. It is also compared with the prediction by Wilson's 

i theory and is found very much improved over the latter. Moreover, the present method is divergence-free 

unlike Wilson's theory. It therefore not only sheds a new light on the nonequilibrium liquid structure of ionic 
solutions and the motion of ion atmosphere in an electric field, but also shows how the Onsager- Wilson 
theory of conductance in electrolyte solutions should be applied to account for the Wien effect on ionic 
conductance. 
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I. INTRODUCTION 



In the preceding companion paper l|, henceforth referred to as paper I or simply I, we have ob- 
tained the exact solutions for the governing equations for nonequilibrium pair distribution functions 
and potentials in a Brownian motion model, and also the solutions for the Stokes equation^- 4], 
which yield the velocity of the countercurrent induced by the motion of ions in a binary strong 
electrolyte solutions under the influence of an external electric field. Therewith we obtained the 
general formulas for the electrophoretic and relaxation time factors that can be made use of to 
define the divergence-free electrophoretic and relaxation time coefficients with which we may cal- 
culate ionic conductance of binary electrolyte solutions. Since the electrophoretic and relaxation 
time factors, f(x, r;£) and g(a;, r;£), obtained are for all allowable values of position (x and r), it 
is necessary to select an appropriate position or a set of positions that must yield divergence-free 
electrophoretic and relaxation time coefficients. Since x and r represent the position of the center 
ion of the ion atmosphere at the given value of external field strength £ we have first determined 
the shape of ion atmosphere from the projection Cp of the axial velocity profile and obtained the 
coordinates of its center position. Then therewith we can calculate the electrohoretic and relax- 
ation time coefficients as functions of external field strength £. The salient point of the procedure 
put forward here is that since the center position (x c , 0) of the ion atmosphere moves along the 
x axis from the origin as £ increases, the electrophoretic coefficient for the center ion must be 
calculated for the moving spherical ion atmosphere at each value of £, which follows the trajectory 
described in Fig. 7 of paper I, and therewith the ionic conductance. In this connection, it should 
be recalled that Wilson [5[ evaluated the electrophoretic coefficient for the center ion located at 
the coordinate origin x = r = for all values of £, which consequently gives rise to a divergent 
term that must be discarded to obtain a finite electrophoretic coefficient. Such a divergent integral 
arises mainly because x = r = is taken for the center ion position for the ion atmosphere even if 
£ > 0. Moreover, his procedure ignores the fact that the ion atmosphere not only is distorted, but 
also migrates in the field direction under the influece of an external field. Therefore the coordinate 
position in the electrophoretic and relaxation time factors must be other than x = r = and also 
a function of £. 

In the present paper, we would like to apply this important point of observation to calculate 
the electrophoretic coefficient and the corresponding relaxation time coefficient, and therewith 
the equivalent ionic conductance in the nonlinear regime of the field strength as a function of 
£; in particular, the Wien effect (Ion equivalent conductance. Thus we would like to show that 
the Onsager- Wilson theory j], u. , is basically capable of predicting in an excellent accuracy 
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a marked deviation from Coulomb's law of conduction if we avoid Wilson's procedure of taking 
x = r = in the electrophoretic and relaxation time coeffficient formulas. When they are calculated 
according to our procedure, the equivalent ionic conductance compares well with experimental 
data on magnesium sulfate solution available in the literature [9]. As a matter of fact, when the 
electrophoretic coefficient f (£) and the relaxation time coefficient g (£) thus calculated according 
to our procedure, the present method becomes not only divergence-free, but also improves the 
accuracy of the Onsager-Wilson theory qualitatively and quantitatively over Wilson's original 
theoretical prediction made with x = r = used for the center position for all values of £. 

The present paper is organized as follows: In Sec. II, the electrophoretic and the relaxation 
time factors are briefly reviewed and the corresponding coefficients are defined. In Sec. Ill, the 
constitutive equations for velocities in terms of diffusion coefficients and mean external forces are 
discussed within the framework of thermodynamics of linear irreversible processes. In Sec. IV, 
the theory of ionic conductance is presented and, upon use of the electrophoretic and relaxation 
time coefficients defined according to the procedure mentioned, the relative ionic conductance is 
calculated as a function of reduced external field strength £ for magnesium sulfate solution in 
excellent agreement with experimental data. Sec. V is for discussion and conclusion. 

II. ELECTROPHORETIC AND RELAXATION TIME FACTORS 
A. Electrophoretic Coefficient 

From the axial velocity v x (x, r; £) completely calculated from the formal Fourier transform 
solutions we have shown the electrophoretic factor f (x, r; £) may be identified: 




(1) 



where f (x, r; £) is given by 
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with £1 (x,r, £), £2 (a^ ©1 c > r j£)> an d ©2 (z^O defined by the integrals 
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It is useful to remind the reader that the second term on the right of Eq. (pQ) is in the same form 
as the electrophoretic effect term — i.e., the last term in Eq. (3) of paper I — in the velocity formula 
obtained on the basis of a heuristic argument. Here the notation is the same as in Paper I, but we 
reiterate that £ = zeX / nkpT is the reduced external field strength and 



Therefore, we see that f (x,r;£) is a measure of variation 1 12] in the external body-force on the 
nonequilibrium ion atmosphere that has resulted from the originally spherical ion atmosphere of 



In Fig. 1, an example of a three-dimensional surface of f (x, r; £) in the case of £ = 3 is plotted 
and the projections of the level curves on the surface are plotted in Fig. 2. The level curves 
appear as quasi-ellipses in the (x,r) plane. In particular, the outermost curve Cp denotes the 
locus of zero of the surface as the intersection with the (x, r) plane. The appearance of two 
quasi-ellipses making up the locus Cp implies that there are two velocity components, axial and 
transversal, developing as the external field is increased from zero. It also means that the spherical 
ion atmosphere formed in the electrolyte solution at equilibrium when the external field is absent 
deforms into two quasi-ellipses, one along the x axis and the other along the r axis transversal to 
the x axis. (This phenomenon in fact gives rise to a Hall effect. It is interesting, but will not be 
a subject of investigation in this paper. We hope to return to this aspect in a future work. Our 
main interest has been and is the quasi-ellipsoid along the axial field direction and its associated 
ion atmosphere in the present series of work.) As is evident from Fig. 1 and Fig. 2, the spherical 
ion atmosphere is now distorted into a quasi-ellipse (a quasi-ellipsoid if rotated about the x axis 
which is the symmetry axis) described by the level curve C p with its center displaced to (x c , 0) 
and tracing the trajectory in the (x c , £) plane shown in Fig. 7 of paper I and also in Table 1, as 
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TABLE I: Trajectory of the center position of the quasi-ellipse as the reduced field strength increases. 



x c 


£ 


x c 


£ 


0.2 


0.5670 


2.4 


0.5960 


0.4 


0.5795 


2.6 


0.5930 


0.6 


0.5900 


2.8 


0.5910 


0.8 


0.5975 


3.0 


0.5895 


1.0 


0.6050 


3.2 


0.5890 


1.2 


0.6070 


3.4 


0.5885 


1.4 


0.6080 


3.6 


0.5880 


1.6 


0.6075 


3.8 


0.5880 


1.8 


0.6065 


4.0 


0.5875 


2.0 


0.6030 


4.2 


0.5875 


2.2 


0.5995 


4.4 


0.5870 



£ increases. Since it would be useful to discuss the feature of f(:r,r;£) in a little more detail, we 
would like to elaborate on it below: 

Since the exact velocity profile consists of a purely mechanical part (v x ) mo and the Brownian 
motion (i.e., stochastic) part (v x ) Brown , which are opposite in their numerical effects in some region 
of space, the electrophoretic factor f(x,r;£), which is a surface in the plane of coordinates x and 
r, has a locus of zero defined by f(x,r;£) = 0; see curve Cp in Fig. 2 of this paper and also Fig. 
6 of paper I. The center of the region bounded by this locus Cp not only migrates along the x 
axis (the field direction) as the field strength increases, but also Cp changes its shape as £ varies. 
When graphically examined, the coordinate x c of the center of Cp traces a certain trajectory as 
£ increases from £ = 0, as shown in Fig. 7 of paper I. Furthermore, f (x, r; £) < if (x, r) is 
within the domain defined by Cp whereas f(x,r;£) > if (x,r) is exterior to Cp. This feature 
implies that the spherically symmetric ion atmosphere at £ = (i.e., equilibrium) with the center 
at the coordinate origin polarizes to a negative velocity domain within Cp and a positive velocity 
domain exterior to Cp in the velocity space while the center of the quasi-elliptic domain of Cp at 
(x c ,0), where x c > 0, is displaced from the coordinate origin where the center of equilibrium ion 
atmosphere was located. Therefore, the spherical ion atmosphere displays two regions of opposite 
signs and the conductance of ions is contributed by ions in the exterior shell, and the electrophoretic 
coefficient should be calculated with the force by the flowing medium acting on the sphere (i.e., 
ion atmosphere) of radius x c , part of which is exterior to Cp, but centered at (x c , 0). This means 
the transversal distance to the x axis (i.e., the radial distance) on the spherical ion atmosphere of 
radius x c (in reduced units used here) in question is r = x c from the coordinate origin (the original 



center). Therefore, the electrophoretic coefficient at £ > for t 



le spherical ion atmosphere at 



(x c ,0) is given by f(£) = f(x c ,r c ;£) = f(x c ,x c ;£), if the Stokes law[2j-|4|] is applied to this spherical 
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ion atmosphere. This, in essence, is the content of the aforementioned procedure for selecting x 
and r in the electrophoretic factor f(x,r;£). This procedure makes the selected values for x = x c 
and r = x c unique because there is only one center and one spherical ion atmosphere of radius 
(l/\/2)^ _1 for each value of £. Note, of course, the set (x c ,r c ) is axially symmetric, because the 
system is axially symmetric around the x axis, which is parallel to the field direction. 

Therefore, the Stokes law should be applied to the spherical ion atmosphere of radius x c in 
reduced units, whose center is located at (x c ,0). Thus, the electrophoretic coefficient at £ is 
defined by 

f(Z) = f(x e ,r e ;S) = f(x e ,x c it). (8) 

It should be noted that f (x, r; £) now is everywhere finite for x > and r > for all £ > 0. We have 
calculated f (£) defined in Eq. ([SD according to the procedure formulated in paper I and mentioned 
earlier. In Table 1, the dependence on £ of x c in the reduced units employed in this work is 
tabulated for the purpose of future use in computing the f (£) and g (£) for binary electrolytes. It is 
interesting to note that this trajectory has a maximum and thereafter diminishes to an asymptote 
as the field strength increases. 

B. Relaxation Time Coefficient 

In paper I, it is shown that the local ionic field e^AX can be obtained from the solutions of the 
governing equations in the form 

ej AX (r) = -^ 9s (i, r; 0) - ^-g (x, r; , (9) 

where the relaxation time factor g (x,r;£) is given by the formula 

(x, r; £) = g c (x, r; f ) + Ag s (x, r; £) (10) 

with g s (x, r;0), g c {x,r;£) and Ag s (x,r; £), respectively, defined by the integrals 

g s (x,r;0) = dye^yh {r^2^p) , (11) 



V2Jo 1 + 2CV J 1 + 2CV 



dy- 



1 f VW+ej e~ xy V 1 + V / l + 2eVUo(wir 



2 7 1 + 2CV 
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It should be noted that Ag s (x, r; £) tends to a constant as £ — )• 0. The integrals in Eqs. (|12p 
and (|13p are also subject to Ineq. (90) or (A20) of paper I. Using this form of relaxation time 
factor g (x, r; £), we have shown in Paper I that Wilson's formula^, 10] g(£) for the relaxation time 
coefficient can be recovered as an approximation if both x and r are taken equal to zero. It is an 
approximation since the contribution from Ag s (x,r;£) must be discarded although its magnitude 
is not negligible when the results of the integrals are evaluated at x = r = 0. 

To be consistent with the electrophoretic coefficient defined earlier, we define the relaxation 
time coefficient g (£) with the relation 

(0 = 9c Oc, r c ; £) + Ag s (x c , r c ; £) = g c (x c , x c ; £) + A# s (x c , x c ; £). (14) 

This relaxation time coefficient a (£ ) will be used in the calculation of ionic conductance as a 

u 

function of the external field strength and, in particular, the Wien effect on the ionic conductance [9, 
fl(j ] of magnesium sulfate solution. In Fig. 3 and Fig. 4, f (£) and g (£) defined in Eqs. ([SD and 
(Q3D are, respectively, calculated and compared with the Wilson's results. Compare these results 
for f (£) and g (£) with the examples for f (x, r; £) and g (x, r; £) for x = r = 0.5 shown in Figs. 8 
and 10, respectively, of Paper I. The values of x and r chosen are arbitrary. 

III. CONSTITUTIVE EQUATION FOR DIFFUSION AND CONDUCTANCE 

Macroscopic description of ionic motions in solution in an external electric field must be sub- 
jected to the laws of irreversible thermodynamics. Since the macroscopic process of interest here 
is steady, we first assume that it obeys steady-state linear irreversible thermodynamic constitutive 
equations — namely, linear thermodynamic force-flux relations. On close examination of the OW 
theory (5] it is evident that the linear irreversible thermodynamic constitutive equations for diffu- 
sion underlie it. If they are found inadequate for constitutive equations in the linear regime then 
we will have to seek a way to extend them to the nonlinear regime. For this, see, for example, 
Refs. 14 and 15. 

to follow a linear constitutive relation 



Here the ionic diffusion flux of ion k is assumec. 
accordin g to the linear irreversible thermodynamics [ill. 



fluids 
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According to the kinetic theory of 



16], for binary mixtures at steady state the thermodynamic force for diffusion is given 



by the linear thermodynamic force 



Xk^Pk ~ — Ffc, (15) 
P 



where x^ = n^/n is the mole fraction, is the partial pressure, p is the pressure, is the mass 
density of species k, and is the external force on the species per mass. In the conductance 
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experiment considered in this work, the partial pressure is uniform in space, so that the pressure 
gradient is equal to zero in first order. In the case of a binary electrolyte solution in which the 
continuum solvent is assumed spatially uniform, there is only one independent diffusion coefficient, 
and the linear thermodynamic force-flux relation or the constitutive equation for diffusion may be 
written as 

Jfc = -uj k pd k , (16) 

where uj^ is the diffusion coefficient of species k. Since the diffusion flux is defined relative to a 
reference velocity by the formula|ll|, 

Jfc = v fc - u, (17) 

where u is the reference velocity suitably chosen, such as the barycentric velocity, which obeys 
the hydrodynamic equations of the solution. There are other possibilities chosen for the reference 
velocity in the literature [ill]. However, in the theory of conductance in electrolyte solutions the 
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161 ] of diffusion is developed, u 



barycentric velocity is generally used, and if the kinetic theory 
is exactly what appears in the equation for the diffusion flux. Since we are generally interested in 
mobilities in the conductance experiments, it is required to calculate v& itself of species k instead 
of the diffusion flux J^. Thus we have 

Vfc = -WfcPd fc + v (r) = 0J k pkF k + v (r) , (18) 

in which we have taken u = v (r), the hydrodynamic velocity calculated previously [1] and presented 
in the previous section. It should be noted that there is no arbitrariness about the identification of 
u with v (r), because generally u in Eq. ()17p is certainly the barycentric velocity obeying the mo- 
mentum balance equation — the Navier-Stokes equation]^ or the Stokes equation^- 4] at the level 
of linear irreversible thermodynamics. By this consideration, we now see that the electrophoretic 
coefficient considered in the theory of conductance is intimately associated with the hydrodynamic 
flow velocity v (r), which is the solution of the Stokes equation for flow of the medium in the binary 
electrolyte solution in the present theory. However, this hydrodynamic velocity to be used is that 
evaluated at (x c ,r c ) at which the electrophoretic coefficient f (x c ,r c ;£) is computed. 

IV. CONDUCTANCE 

The velocity of ion k at (x c , r c ) is computed relative to a suitable reference velocity — for 
example, the solution of the Stokes equation. Identifying pk^k with e^Xf = efcX+e^AX and taking 
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the reference velocity u(x c ,r c ) with the axial velocity v x at (x c ,r c ), we find 

ft\ v ^kekV'K . . \z k \eXK . 
Vfc (0 = u k e k X __g (f) - _^_f (f ) (k = i,j). (19) 

Here we have taken the ionic field AX as for the relaxation time effect in Eq. ([SJ 

ejAX = -^g (xe, r c ; £) = -^g (0 (20) 

without the field-independent term which has nothing to do with the relaxation of distorted ion 
atmosphere, because the relaxation time coefficient is representative of the relaxation of asymmetry 
of the ionic atmosphere. Similarly, the electrophoretic effect is represented by 

/ ^\ Z€-X K „ . j_ x ZcX K . . . . 

v * ® = -7T7*—f ^ rc! $ = -T7^— f (?) (21) 
by Zirr/o by Ztttiq 

without the field-independent term. 

The mobility in electrostatic units is then given bvfic]] 

„ - Vfc (^) / 99 A 

u k - x , (22) 

and it may be written in practical units as 

1 / u) k e k [i'K \e k \K \ 

u k = — [\e k \ u k - -^-0(0 - ^f(0j • (23) 

At infinite dilution, k — >• and the limiting mobility is given by 

u° k = i^*. (24) 

fe 300 v 7 

Since the limiting equivalent conductance is 

A° fc = Ful (25) 

where -F is the Faraday constant (F = 96493.1), we have the equivalent conductance of ion k in 
practical units given in the form 



ze 



\2 



kA9 F|efe|/« 



A ' = A! -w ,,!| -^» ft (26) 

Define the equivalent conductance of the binary electrolyte 

A = A + + A_ =A!+A 2 . (27) 



Then the equivalent conductance for a binary electrolyte is given by 

(zefnA _ F(\ ei \ + \e 2 \), 
2Dk B T m> 6 x 300^27^0 



A (0 = A° - %^r0(O - IM+Ml m (28) 
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with A = A° + Ag. 

In experiment, the relative equivalent conductance 

AA(C) = A(O-A(0) (29) 

is reported as a function of the field strength £. Since in the limit of £ = 

we obtain the relative equivalent conductance 

This formula is computed for various values of reduced field strength £, and AA(£) /A (0) 
calculated is plotted against £ and compared with experimental data[s, 10] on a MgS04 solution in 
Fig. 5. It is the only data of this kind available in the literature for strong binary electrolytes which 
we can use for comparison on the Wien effect. We have taken d = 2i?MgS0 4 = 2 x 3.12 = 6.24 for the 
diameter of MgS0 4 , T = 291 K, D = 81, c = 1.22 x 10~ 3 mole/liter, r] = 0.0105 poise, the limiting 
conductance values flo| Aq = 53.06 for Mg +2 and Aq = 80.0 for SOJ 2 , respectively, and 120.36 
g/mole for the molecular weight for MgS04. In Fig. 5, the solid curve is the prediction by the 
present theory, the solid squares are experimental data, and the broken curve is the prediction by 
Wilson's method, that is, the prediction with the electrophoretic coefficient /(£) and the relaxation 

. Note that they are tabulated 

in the monograph of Harned and OwenjlQ]. The OW theory prediction is much smaller than the 
experimental conductance. Furthermore, the low field (£ < 1) behavior predicted by the Wilson's 
result appears to be markedly different from the experimental data and the present theoretical 
prediction, the latter being almost linear with respect to £, while Wilson's prediction appears to be 
nonlinear, i.e., £ a with a > 1. We believe that the procedure used for selecting the x and r used for 
calculation of the electrophoretic and relaxation time coefficients and thus the ionic conductance 
appears to be well supported by the comparison made. 



V. DISCUSSION AND CONCLUDING REMARKS 



Ionic liquids are ubiquitous and it is important to understand their physical behaviors in phys- 
ical, chemical, and biological systems under the influence of external electromagnetic fields. On- 
sager's theory|7] of conductance is an important benchmark in the development of the theory of 
transport processes in ionic solutions, although his theory appears to have been almost forgotten 
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17]. We believe it is unwarranted 



about and paid little attention these days except one recent work 
to ignore his basic ideas on the subject matter because the theory seems to provide much valuable 
insight for further development of theories of transport processes in ionic solutions and plasmas. 
In the present series of papers, we have revisited and examined the details of his theories and 
have been able to obtain full evaluations of the formal results of the OW theory. The preceding [l| 
and present papers describe our understanding of the inner workings of his theories and their ap- 
plication to an experimental system. We have found that the solutions of the OF equations[18], 
the Poisson equations [3], and Stokes equation^ for the velocity of the medium would require a 
more mathematically complete treatment than originally implemented in Wilson's dissertation^, 
especially, to avoid divergence difficulty that was unfortunately not dealt with squarely. In paper 
I[l] we have exactly evaluated the formulas for the axial and transversal velocities in a binary elec- 
trolyte solution, which may be deemed one of the goals of the theory of ionic conductance. In the 
present paper, we have applied the axial velocity formula to calculate the equivalent conductance of 
a binary strong electrolyte, having elucidated the appropriate set of (x, r) to use by a computation- 
based procedure to select a unique set of (x, r) as a function of £. Although this procedure lacks 
a precise theoretical foundation based on irreversible thermodynamics or hydrodynamics, it does 
provide us with an unambiguous procedure to obtain correct estimates of conductance over the 
entire range of £ experimentally studied. In this sense, we may conclude that it is a practical and 
useful procedure for understanding the experimental data. In any case, the procedure is founded 
on the observation of how the ion atmosphere behaves in the applied external field. The study of 
its behavior with respect to the field not only provides much insight into the interplay between the 
long-range Coulombic interactions between charged particles and the interaction of the ion atmo- 
sphere with the external field producing a rather intricate ionic liquid structure, but also shows 
how charged species diffuse in such a liquid under the influence of the applied electric field. This 
picture of the ionic liquid behavior that we have gained in the present paper owes to the fact that 
the Brownian movement model, namely, the solutions of the governing equations, provides us with 
exact nonequilibrium distribution functions, ionic potentials, and velocity profiles that can be eval- 
uated rigorously either numerically of analytically. This latter feature is an indispensable aspect 
of the OW theory of ionic solutions we have shown possible in this work. We emphasize that the 
present work disposes of the divergence difficulty inherent to Wilson's procedure [5i], which arose 
from an inadequate treatment of the Fourier transform integrals for the solutions, and the theory 
of Wien effect is thereby made practically amenable to theoretical investigation and satisfactorily 
accountable theoretically. 
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The theoretical significance of the Onsager theory of electrolytic conductance [7;] lies in its abil- 
ity to take into account the non-linear field dependence in the form of local body-force that is 



obtained through the Onsager-Puoss equations 



son equations 



181 ] for the pair distribution functions and the Pois- 



19[ | which are coupled together. The local body-force thus obtained is then made 



use of in the Stokes equation for the velocity of the medium subjected to the external electric field. 
The present paper enables us to conclude that the Onsager theory is now made capable of treating 
conduction phenomena of dilute electrolyte solutions irrespective of the applied field strength and 
predicts a correct behavior qualitatively and quantitatively in contrast to the conclusion drawn in 
the literature in the past that it predicts too small values for equivalent conductance [lo| . Thus 
it, in fact, appears to provide a great deal of insights into diverse electrical conduction phenom- 
ena in subjects other than traditional electrolyte solutions studied in physical chemistry, which 



are lately studied under different guises and subject names, such as semiconductors 



liquids 



22, 



231 ]. ion channels 



plasmas 



20 



211 ] . ionic 



251 ]. etc., but still involve charged species subjected to an 



external electric field gradient, often rather strong because of small system sizes commonly stud- 
ied in recent years. Therefore the basic physics underlying the phenomena in the aforementioned 
subject fields is unchanged other than sizes of the systems and conditions, and in this respect 
Onsager's theory can give us a valuable lesson on how to treat ionic solutions and ionized gases in 
external electromagnetic fields. 

In this work and also in the OW theory Q], the treatment of the thermodynamic force-flux 
relation is still that of linear irreversible processes, that is, the constitutive relations for ionic motion 
and flow of the medium are those of the theory of linear irreversible processes ll| . In addition to 
this restriction, it is a theory limited to sufficiently dilute electrolyte solutions. Therefore these two 
important restrictions must be removed for a more comprehensive theory of electric conductivity 
in experimental systems of current interest. We leave these aspects to future studies, but just refer 
to reader to Refs. 14 and 15 for nonlinear constitutive equations. 
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Figure Captions 

Fig. 1 The electrophoretic factor f (x,r, £) is plotted in the (x,r) plane in 3D in the case of 
£ = 3. 

Fig. 2 The level curves of the surface f(x,r,£) in Fig. 1 is projected onto the (x,r) plane. 
The level curves form two sets of quasi-ellipses with the major axes on the x and r axes, respec- 
tively. The outermost level curve Cp represents the zero of f (x,r, £), namely, f (x,r, £) = 0. This 
outermost level curve Cp indicates how the spherical ion atmosphere assumed by the equilibrium 
ion atmosphere when £ = migrates from (x, r) = (0, 0) to a point on the positive x axis and also 
is distorted to a quasi-ellipse with the major axis lying on the x axis by the action of the external 
electric field as £ increases, and similarly to the quasi-ellipse with the major axis lying on the r 
axis perpendicular to the x axis. The center of the curve Cp follows a trajectory described by the 
curve shown in Fig. 7 in paper I. The trajectory grows from zero at £ = to a maximum and then 
decreases to a plateau value as £ increases. See also Table 1 of the present paper. 

Fig. 3 Comparison of the electrophoretic coefficients of the present and Wilson's theories. The 
solid line: the present theory; the dotted line: Wilson's result. 

Fig. 4. Comparison of the relaxation time coefficients of the present and Wilson's theory. The 
solid line: the present theory; the dotted line: Wilson's result. 

Fig. 5. Relative equivalent conductance compared with experiment and the prediction by 
the Onsager-Wilson formula. The solid symbols are the experimental data on MgS04 solution, 
the solid curve the present theory, and the broken curve the prediction by the Onsager-Wilson 
formula. It should be noted that the behaviors of the relative conductance of the present and 
Wilson's theories are qualitatively and also quantitatively different. 
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